ON THE DIRICHLET PROBLEM FOR FIRST ORDER LINEAR HYPERBOLIC 
PDES ON BOUNDED DOMAINS WITH MERE INFLOW BOUNDARY 
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Abstract. Here we study the Dirichlet problem for first order linear and quasi-linear hyperbolic PDEs 
on a simply connected bounded domain of R 2 , where the domain has an interior outflow set and a mere 
inflow boundary. By means of a Lyapunov function we show the existence of a unique solution in the 
space of functions of bounded variation and its continuous dependence on all the data of the linear prob- 
lem. Finally, we conclude the existence of a solution to the quasi-linear case by utilizing the Schauder 
fixed point theorem. This type of problems considered here appears in applications such as transport 
based image inpainting. 
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1. Introduction. The subject of this paper is hyperbolic linear and quasi-linear 
scalar PDEs of first order in two space dimensions. We consider the Dirichlet prob- 
lem on simply connected bounded domains O and bounded functions uq of bounded 
variation as boundary data. The linear PDE is stated in the space BV - the functions 
of bounded variation - and formulated as 

(c(x),Du) =f{x)-C 2 in n\S, u\ dn = u . (1.1) 

Here, ( . , . ) denotes the scalar product in R 2 , Du is the derivative measure of u 6 BV 
and £ 2 the Lebeseue measure on R 2 . 

Because the PDE (1.1) is hyperbolic and our aim is to ensure global existence and 
uniqueness, we have to rule out the case of characteristic points and therefor need 
an additional requirement. Our additional requirement is the existence of a global 
Lyapunov function T : Cl — » R for the transport field c. That is c is assumed to satisfy 
the Lyapunov or causality condition 

(c{x),VT(x)) > B- |VT(x)| , B>0. (1.2) 

Besides that, the Lyapunov function T is zero on the boundary and increases towards 
its maximal level E. 

Since the boundary is a level line of T, condition (1.2) means in particular that 
every point of the boundary is an inflow point. On the other hand PDE (1.1) is an 
advection equation, thus we need an outflow set. For our case the outflow set is in 
the interior of the domain and is given by L. The outflow sets considered here are 
connected tree-like sets compactly contained in O. 

Assuming c G C 1 (n\E) 2 , equation (1.1) is a well-behaved advection equation 
and thus we will employ the method of characteristics to construct the solution. A 
related problem, the Cauchy problem for quasi-linear equations in conservative form 

dtu(t,x) + d x a(t,x,u(t, x)) = , u(0,x) = uq(x) , (1.3) 

with bounded data uq of locally of bounded variation, has already been studied, see 
e.g. [ ] and [ ]. Here, the non-linearity of a(t,x,u) w.r.t. the variable u causes the 
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solution to develop shocks within finite time, which in turn makes a notion of gen- 
eralized solution necessary. In [ ], the author proves the existence and uniqueness, 
under some suitable assumptions on the jump discontinuities, in BV D L°° employing 
the method of vanishing viscosity. 

Our problem is slightly different. The set E is globally attractive by the causality 
condition (1.2) and thus characteristics will meet there. Because equation (1.1) is a 
well-behaved advection, the solution cannot develop shocks before the characteris- 
tics reach the set E but we will in general obtain shocks on E. This is the case even 
for smooth boundary data. 

In this paper we will show existence, uniqueness and stability of the solution to 
problem (1.1) in the space BV(Cl) n L°°(fi). 

An important feature of our solution is that it is defined on all of O. That means, 
we obtain a description of the shocks at E in the B V-framework and so, in this sense, 
we can close the gap: u solves problem (1.1) in D\E, but u G BV(fi) and not only 
u e BV(Q\E). 

The benefit of "closing the gap" is the good functional analytical properties of 
BV(d), in particular compactness, which we will exploit to tackle the quasi-linear 
case: 

(c[u](x),Du) = f[u](x)-C 2 in Q\E, "Ian = "o ■ (1-4) 

In contrast to equation (1.3), we consider quasi-linear problems in non-conservative 
form but we allow the dependence on u of the transport field c and the right-hand 
side / to be of a general functional type. This means that they depend not only on 
the value u(x) but on the whole function, thus the coefficients of the PDE (1.4) are 
maps 

f:X^Q lr u^f[u], /[.](*) :X^R , u ^ f[u](x) , 

9 (1-5) 

c : X — > Qi , u — y c[u] , c[ . } (x) : X — > ]R , u — > c[u] (x) , 

with X C BV(d) and Q\, Qi being subsets of suitable function spaces defined on 

n\E. 

We will approach the quasi-linear problem (1.4) by using the theory of the linear 
problem. Fixing the functional argument of the coefficients by some v £ X, we obtain 
the linear PDE 

(c[v](x),Du) =f[v](x)-C 2 in Q\E , w|an = "o, (1-6) 

with a solution U[v] EX depending on v. If now the operator U admits a fixed point 
u — U[u], we obtain a solution to problem (1.4). 

The compactness mentioned above and the stability properties of the linear case 
allow us to establish the existence of a solution u = U [u] by employing the Schauder 
fixed point theorem. 

An application of the quasi-linear equation (1.4) is transport based image in- 
painting. The term image inpainting means the retouching of undesired or damaged 
portions of a digital image. Telea suggested in [ ] a fast algorithm for digital image 
inpainting which is however not adapted to the image under consideration. 

In [ ], we have analyzed Telea's algorithm in a continuous setting and have seen 
that it is consistent with the following model 



(VT(;t),Vu) =0 in n\E, 



u \dd = u 
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Moreover, we have generalized this linear model to problem (1.1) with side condition 
(1.2) and finally we adapted the coefficients of the PDE to the image which resulted 
in a quasi-linear model of type (1.4). With our image adaption we obtained an im- 
age inpainting method which produces results of high visual quality and performs 
almost as fast as Telea's algorithm. 

In [ ], we left open the question for the well-posedness of our inpainting model. 
This question is positively answered in [6] as an application of the theory presented 
in this paper. 

Outline of the Paper. Section 2 is devoted to the full description of the linear prob- 
lem (1.1) and its requirements. The considered outflow sets E are connected tree-like 
sets compactly contained in (1, thus the biggest part of the requirements is concerned 
with the complicated structure of E in order to define Lyapunov functions T and 
transport fields c which are well-behaved close to and on the outflow set. In sec- 
tion 3 we discuss the non-linear coordinate system induced by the characteristics, 
in particular the features which are crucial for integral transformations later on. In 
section 4, we construct a solution u £ BV (fi) by using the method of characteristics 
and derive a priori bounds on the norm || u \\ gv(n) which depend only the data of the 
problem and features of the Lyapunov function T. The main result is theorem 4.1. 
Section 5 is concerned with the uniqueness and the continuous dependence of u on 
the data. Special attention is given to the non-linear dependence of the solution on 
the transport field c and on the Lyapunov function T. The main result is theorem 5.1. 
In section 6 we turn to the quasi-linear problem (1.4). Exploiting the theory of the 
linear problem, in particular the continuous dependence on the data and that we can 
close the gap as u 6 BV(fi), we conclude the existence of a solution by the Schauder 
fixed point theorem. 

2. The Linear Problem and its Requirements. The purpose of this section is to 
summarize all the requirements on the data of the linear problem: 

(c(i),D«}=/(x)-£ 2 in n\E, 

"|an = "0/ (2-1) 
(c(x),VT(x)} > 6- |VT(x)| . 

The first set of requirements is concerned with the domain O, the outflow set E and 
the Lyapunov function T. Later, in section 4, we will construct the solution by the 
method of characteristics. For the characteristics we will see that their time variable 
can be identified with the values of the Lyapunov function T and that the outflow set 
E is the location where characteristics end. For these reasons, we call the Lyapunov 
function T time function and the outflow set E stop set. The second set of require- 
ments is concerned with the transport field c, the right-hand side / and the boundary 
data uq. 

Requirement 2.1. (Domains) Domains O c IR 2 are assumed to be open, bounded 
and simply connected and to have C 1 boundary. 

Because of requirement 2.1 the boundary 3Q of a domain is a simple closed C 1 
curve. By 7 : ]R — » dCl we denote a generic regular and periodic parametrization of 
dd. That means 7 G C 1 (1R, dD.) is surjective and 7'(s) 7^ Vs £ IR. Furthermore, by 
I = [a, b[ C JR we denote an interval such that 7)7 is a generator of 7. 

Next, we collect the properties of admissible time functions. For our problem 
time functions are global Lyapunov functions whose range corresponds to a finite 
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time interval. As time is usually a positive quantity which increases, we define the 
stop set to be the maximal level of T while in literature the stop set is often the 
minimal level, see e.g. [1]. 

Requirement 2.2. (Time functions) Time functions are of type T : CI — >■ K. The 
upper level-sets of T are denoted by 

Xt>a := {x e n : T(x) > A} . 

We assume that time functions T satisfy the following conditions: 

1. T G c(n). 

2. Tfte boundary off! is the start level: T\^ =0. 

3. T incorporates a stop set E iw'f/z stop time 1; 
a) T(x) < 1 x 6 Q\E. 

W T|s = 1, z.e., E zs zTze maximal level ofT. 

4. T increases strictly from dD, towards E. TTzaf means that any upper level-set Xt>a zs 
simply connected and 

Xt>a=X¥^ VAe[0,l[. 

5. Any proper upper level-set is a future domain: for every A e [0, 1[ the set Xt>a satisfies 
requirement 2.1. Furthermore, the field of interior unit normals to the \-levels 

XT=A=dXT>A, Ae[0,l[ 

ofT is denoted by N : D\E — > S 1 . N is required to be continuously differ entiable and 
extendable onto dCl, i.e., N 6 C 1 (C\E). 

6* T 6 C 2 (n), icxt/z VT(x) =0»ieL. 

Remark on 6.*: this assumption is in order to ease things in the passages that 
follow. Because of 6.*, we obtain a simple description of the field N on fi\E by 
N(x) = VT(x)/|VT(x)| and N is continuously differentiable and extendable onto 

an. 

In part 3 of requirement 2.2 we have assumed that T features a stop set E. Here 
we state the geometric properties of allowed stop sets. 

Requirement 2.3. (Stop sets) We assume that the stop sets E satisfy the following 
conditions: 

1. His a closed subset of CI. 

2. E is either an isolated point, or a connected set with tree-like structure (no loops). 

3. If E is not an isolated point, it consists of finitely many rectifiable C 1 arcs Ej-: 

n 

E = (J Z k . 

k=l 

The collection {Hk}k=i,...,n !S assumed to be minimal in the number n of arcs, so E is 
decomposed by breaking it up at corners and branching points. 

Furthermore, we require for each arc Ejt that its relative interior Ej. has a given orientation 
by a continuous unit normal : Ej- — >■ S 1 . 
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In the case in which E is not just an isolated point, we also need good behavior 
of the maps T and N close to and on the stop set E. For this purpose, we use the 
following concept of one-sided limits towards z E Ej.: 

If x E B r (z)\E and if r > is sufficiently small, the projection p of the point x is 
unique. In view of this feature we say a point x E B r (z)\L is on the right-hand side 
or plus side (respectively, on the left-hand or minus side) of E^ if 

i^T = +«*(P) ( T^T = -MP)) ■ (2.2) 



Therewith, a sequence (x„)„ GN , x n E Q\E, tends to z 6 Ej. coming from the plus 
side (respectively, the minus side), in symbols 

x n -> z + (x„ -> z_) , (2.3) 

if the sequence converges towards z and almost all elements x„ are on the plus side 
(respectively, minus side). 

With the concept of one-sided limits, the good behavior of the maps T and N is 
summarized in the following requirement. 

Requirement 2.4. (Good behavior at E) 

1. Requirements on T: 

Let y E E and h E S 1 . Let p = p(y,h) be the best possible order for the asymptotic 
formula 

T(y + rh) = l-0(rP) , r -> 0+ . 

We require that there is a bound q such that sup sup p{y,h) < q. 

yes |h|=l 

2. Requirements on N: 

a) N has one-sided extensions onto the relatively open components Ej. and those exten- 
sions are given by ±.n^: 

N+(y) := lim N(x) , N+(y) = -n k (y) , 

N-{y) : Jim \';v:, N~ (y) = n k {y) 

for every y E Ej-. 

b) The derivative DN has one-sided extensions onto the relatively open components E fc , 
i.e., 

(DN) + (y) := _Km DN(i) , (DN)~(y) := _Um DN(x) 

exist for every y E Ej.. 

c) |DN| 6 L 1 (O), i.e., poZes of |DN| at corner-, branching- and terminal nodes o/E are 
integrable. This feature is assumed to hold in the case in which E is an isolated point 
as well. 

So far we have the requirements on domains, stop sets, and time functions. Now 
we turn to the assumptions on transport fields. Here, for a given time function T 
with stop set E, the causality of the transport field w.r.t. T and its good behavior 
close to E are the main concern. 
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Requirement 2.5. (Transport fields) Transport fields of type c : Q\£ — > M 2 are 
required to satisfy: 

1. c s C 1 {Cl\L) 2 and c features the following properties: 

a) c and Dc are continuously extendable onto dfl. 

b) IfL is not just an isolated point, then c and Dc have one-sided limits on the relatively 
open C 1 arcs £ fc ofL: 

c + (y) = lim c(x) and c~(y) = lim c(x) , 
(Dc)+(y) = lim Dc(i) and (Dc)"(y) = lim Dc(z) , 

for every y G E*. 

2. Ura'f speed and causality condition: 
a) \c\ = 1. 

W TTiere is a lower bound B > suc/i f/zaf 

P < (c(x),N(x)} < l Vien\L. 

cj Conditions a) and b) hold for the one-sided limits as well, i.e., 

\c + (y)\ = \c-(y)\=l 

and 

B<(c + (y),N + (y)) <1 , p<< C -(y),N-(y)><l, 

whenever y belongs to some Lj.. 

3. Let Zfc, k 6 {1, . . ., m} denote the terminal-, branching- and kink nodes of E. For every 
e > such that each disk B £ (zj c ) is compactly contained in O, we define the set 



V £ :=E U \jB £ (z k ). 
k=l 

a) For every admissible £ > 0, there is a bound M £ such that 

\Dc{x)\<m £ , Vxen\y £ 

b) \Dc\ eL 1 (Cl),i.e.,polesof\Dc\atz k ,ke{l,...,n}areintegrable. 
Finally we need a right-hand side and boundary data. 

Requirement 2.6. (Right-hand sides and boundary data) We assume right-hand sides 
f : fl — > 1R and Dirichlet boundary data uq : 3D — > Ktobef 6 C 1 (n)and«o £ BV(9Q). 



Remark: with «o G BV(9n) we mean that for every regular periodic parametri- 
zation 7 of dd, the pull-back 7*»o = u 7 is a periodic BV function. Moreover, 
because dD, is one-dimensional, 7*«o is BV function of one variable. Therefor the 
boundary data is essentially bounded. 

In the following sections we assume, if not explicitly stated otherwise, that do- 
mains O, stop sets E, time functions T, transport fields c, right-hand sides /, and 
boundary data uq satisfy the requirements above. 
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3. A Customized Coordinate System. The requirements of the previous section 
in combination are such that the family of characteristics which is the solution of the 
initial value problem 

y> = c(y) , y(0) = x G fi\E , 

gives us a customized coordinate system for problem (2.1). The subject of this section 
is the features of this non-linear coordinate system. The first lemma is about the 
boundedness of the arc-length of characteristic curves. This property is important 
because the ||.|| By-norm bound on the solution of (2.1) depends on it. 

Lemma 3.1. 

i 

a) Let q be the bound from requirement 2.4 part 1 , let cp(t) := —ti and let 

T (x) :=l + <p(l-T(x)) . (3.1) 

Then the gradient VTo of the transformed time function blows up at E and is bounded 
below 

|VT (x)| > m > Vi6fi\S. 



b) For every regular C 1 curve x : [0,a[ — > Q\E (a = oo admissible) that satisfies the 
condition 



0</J<^M*(T)),p^), Vtg [(),«[, 
the arc length of x is uniformly bounded by 



(3.2) 



1 

arclength(x) < . (3.3) 

j6-m 

Proof. 

a) The function To is well-defined since < T < 1 and its derivative is 
VTo(x) = - T(x)) ■ (-VT(x)) =: H(x) ■ , 

with 

H{x) := -(1 - T(x)) V . |VT(x)| > , xe n\S . 

Let y 6 E, h G S 1 and r > 0. Then, by requirement 2.4 part 1 we have 
l-T{y + rh)=C 1 rP, p = p{y,h) , Q > . 

Because T G C 2 (D) and VT|^ = we have 

\VT(y + rh)\=C 2 rf- 1 , 
with the same p as before. Combining these asymptotic formulae, we obtain 

H(y + rh) = C 3 r i r^ 1 = C 3 r i . 
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Since, by requirement 2.4 part 1, q > p{y,h) holds uniformly, we get a blow-up 
lim |VT (y +rh)\ = lim H(y + rh) = oo , 

r—^0 + r->0+ 

for any choice of h £ S 1 and every y £ E. 

We show next that | VTo| > mo > 0. Assume by contradiction that 
inf H(x) = inf |VT (x)| =0. 

and choose an open neighborhood U of E such that H\ u > M, for some constant 
M > 0, which is possible because of the blow-up. Then the restriction onto the 
compact complement H = H\^ u , being a continuous function, must take the 
minimum 

min H(x) = 

x€Ci\U 

at some point x G Ci\U. But then, the definition of H implies 

H(x)=H(x)=0 ^ |VT(*)|=0, 

which is a contradiction, since $ £ E. Thus, H = | VTo| has a minimum wzo which 
is greater than zero, 
b) Using Tq we estimate the arc length from above by 

f 

1 > 7b(jc(f)) - T (x(0)) = j (VT (x(t)),x'(t)) rfr 



= / (N(x(r)),^^y\WT (x(r))\-\x'(r)\dr 

t t 
>pj |VT (x(T))|-|x'(T)|dT > p-moj \x'(r)\dr. 


The limit t — »■ a finally yields the uniform bound on the arc length of such curves 

1 

arclength(x) < j— , 
which depends only on f> and information from T. 

□ 

Because of its nice properties, the transformed version To defined by equation 
(3.1) will be identified - instead of T - with the time variable of the characteristics. 
This is the subject of the next lemma. Whenever we speak about Tq we mean this 
transformed version of a given time function T. 

Lemma 3.2. 
a) The initial value problem 

y' = c(y) , y(0) = x e Q\E , (3.4) 
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has a unique maximally continued solution y : ]f_, t + [ — > IR 2 , with -co < f_ < < 
t + < oo. 

Every trajectory y connects the sets dCl and E, i.e., 

lim 1/(0 e , lim y(£) e E ■ 

For ez?ery pomf z e ^ m f/ze relative interior of some C 1 arc o/E, there are exactly two 
trajectories which hit z in the limit t — > t + , one for each side of Ej.. 

b) The solution yo of the velocity-transformed IVP 

c 

j/ = c (y), y(0)=xeO\E, c := ^— , (3.5) 

satisfies T (yo(£)) = i + Tb(x). 
Proof. 

a) We consider the IVP (3.4). Because c is Lipschitz continuous by requirement 2.5 
part 3a) there exists a maximally continued, unique solution y with time domain 
](_,£+[ and € ]t-,t+[. 

Because of the unit speed condition \c\ — 1, y never stops inside Q\E and never 
blows up. The causality condition (requirement 2.5 part 2b) ) implies, by 

^T (y(0) = (Vr (y(f)),c(y(f))> > m ■ f> > , (3.6) 

that To(y(f)) strictly increases at least at a rate of m^ ■ fi. 

Thus, y collapses at the boundary of Q\E and by (3.6) there follows: 

• Going forward t — »■ t+: collapse at E after finite time f+ , 

• Going backward f — > f_: collapse at after finite time f_. 

Because of unit speed the values t+ , f- are exactly the arc lengths, which are 
finite by lemma 3.1. 

Assume now that z £ and consider the "plus side", i.e., the side where n^{z) 
points. Since c and Dc both extend from the plus side onto £^ by c + and (Dc) + , 
the backward IVP 

y' = - c {y), y(0)= z > with c(z) := c+(z) . 

has a unique solution that starts atz 6 and evolves away from E into the 
plus side. Hence, vice versa there is only one solution y of the forward IVP (3.4) 
that comes from the plus side, heads for z 6 E/ c , and hits z in the end. The same 
argument holds for the minus side. 

b) We consider again the forward IVP (3.4), the initial value X of which satisfies 

Tq(x) < 1 — maxTo(z). For A with Tq(x) < A < 1, there is a unique time t(A), 

zefi 

when the solution y of (3.4) crosses the A-level of To. This is because To(y(f)) 

increases strictly by (3.6), so y crosses the A-level only once . 

Then, viewing T as a function of A, the implicit function theorem applied to 

r (y(T)) = a 

yields the differentiability of T w.r.t. A and the derivative: 

1 



t'(A) 



(VT (z),c(z)) 



!=y(r(A)) 
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Using T, we now change the independent variable 

y (A) :=y(T(A + A )) , 

where Ao := T}(x) and t(Ao) = 0. Then, yo satisfies the initial condition yo(0) = x 
and has the derivative 



y (A) = y'(x(A + Ao)) • r'(A + A ) = (c(z) ■ Jyj^^ ) 

co(yo(A)) • 



z=y(r(A+Ao)) 



(VT (z), C (z)) 



z=y (A) 



Consequently, yo is the unique solution of the transformed IVP (3.5) and — by 
construction — satisfies 

T (y (A)) = T (y(x(A + Ao))) = A + A = A + T (x). 

□ 

By lemma 3.2 part b), we get the useful property that - when using cq instead of 
the original transport field c - the time variable of a characteristic yo is given by Tq. 
We use this feature to introduce new coordinates on f]\L whose concept is similar 
to polar coordinates on a disk: the level lines of To will play the role of the concentric 
circles and the characteristic curves that of the radial lines. 

Whenever we speak about cq we mean the transformed version of a given trans- 
port field c according to lemma 3.2 part b). 

COROLLARY 3.3. (Polar coordinates) Let jbea regular periodic parametrization ofdd 
with generator y\ i, I = [a,b[. 

a) The general solution f (f,s) of the forward IVP 

y'=c {y), y(0)=7(s), 
defines a diffeomorphism £ : ] 0, 1 [ X ] a, b [ — > Q\ (L U S) , where 

S = {£{t,a):t £]0,1Q . 

b) Let rj(t, x) denote the general solution of the backward IVP 

y' = -c (y), y(0) = xen\(s:us). 

Then the inverse map (x) = (t(x),s(x)) T is given by 

r 1 (x) = (T (x), 7 - 1 (^(T (x),x)) ) T , ien\pus). 
The relation between £ and n is 

g(t f s(x))=ti(T (x)-t,x). (3.7) 

c) If 7 is oriented clockwise, the Jacobian D£ = (3f£|9 s £) of the diffeomorphism £ has a 
positive determinant and the estimate 

0<detD£<|9 t £||a s £|<^^ (3.8) 

holds true. Ify is counter-clockwise, the assertions hold for — detD£. 
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d) For each of the relatively open C 1 arcs E^ of E, we can find - w.r.t. the orientation n^ - 
two subsets fa + and fa of I such that the maps £(l,s) with s e fa + and £(l,s) with 
s 6 fa _ are froffa regular C^-parametrizations of £/ c . 

Proof. The proof of these statements can be found in [6, section 3.2]. □ 

Our customized coordinate system is now given by the diffeomorphism £ of 
corollary 3.3 part a) and we can turn to solving problem (2.1). 

4. Existence of a Solution. In this section we construct a candidate solution to 
the linear problem (2.1) and study its properties. The tool for construction is the 
method of characteristics (see e.g. [5] or [ 1 0]) which we apply it to the scaled PDE 

(c (x),Du) =f (x)-C 2 in n\E, (4.1) 

with 

cix) fix) 

co(x) := , , s ' , xx and f (x) := . , . J „L , ., . (4.2) 

For PDE (4.1) the characteristic equation is exactly the forward IVP from corollary 
3.3 a) 

y' = c {y), y(0) = 7 ( s ) / 

and the family of forward characteristics is £. Clearly the solution n of the corre- 
sponding backward IVP is the family of backward characteristics. 

Let now v(t, s):=uo£(f, s). Then - at least formally - the partial derivative oiv 
w.r.t. t is given by 

d t v(t,s) = (Vwo£(f,s),af£(U)) = (VM,c )o^(f,s) = / O o£(f,s), 

together with the initial condition 

v(0,s) = u(Z(0,s))=u(y(s)) = 7 *Uo{s) . 

Here, 7* denotes the pull-back operation. Thus, by the fundamental theorem of 
calculus, we obtain 

f 

v(t,s) = 7 *u (s) + J / o£(T,s)dT =:v l (t,s) + v 2 (t,s) . (4.3) 


The function v represents our candidate solution u in characteristic variables (f,s). 
By using the inverse map from corollary 3.3 b) and the relation (3.7) between £ 
and rj, we push v forward onto Q\E to have u = v o in the original variables x, 

T {x) 

u(x) = u (t](Tq(x),x)) + y foOT](T,x)dr -. u\{x) + u z (x) . (4.4) 

o 

Equations (4.3) and (4.4) present our candidate solution. The next theorem shows 
that the candidate belongs to the space BV(fi) n L co (n). 

THEOREM 4.1. (Element ofBV) The candidate solution ufrom (4.4) with its decompo- 
sition u = u\ + «2 has the properties: 



12 



T.MARZ 



a) u is an element ofBV(Cl\L) n L°°(n). Its L°°(0) norm is bounded by 

Hull < llu II + ll/l|o ° 

ll"llL~(n) < l|MollL-(an) + j-^ ■ 

The derivative measure ofu is 

Du — Cq(x) ■ ]i + X7u2(x) ■ C 2 with u := £j ^C 1 <8> Dy*uoj , 

whereas the total variation is bounded by 

m <_ Mnv :- !^ + (IZb + j^).< ?( n, 

+ ^-( l|D<r|| y<™ + l|DN|1 ''™) ' 

b) u extends onto E, i.e., u is an element ofBV{Cl) n L°°(Q). In comparison to part a) the 
extension introduces - in the derivative Du-an additional jump part for every C 1 arc 5^ 
ofZ: 

n 

Du = cj-(x) • + Vu 2 (x) ■ C 2 + YK u tS x ) ~ M z t ( x )) -T^LZk , 



'-k ' 

k=l 



where w Ejt and are Ze/it and rigM interior BV -traces ofu on H k . The bound on the 
total variation is added up by 

\Du\(CL) <M n \ E +2- ||«|| L - (n) -H 1 ^) . 



Proof. Before going into the details of the proof, we summarize some facts con- 
cerning the change of variables. If <p 6 C] (Q) or cp 6 C\ (0\E) is a test function, we 
will denote by ip(t, s) := <p o £(£, s ) the test function in characteristic variables. By the 
chain rule we then obtain the derivative with respect to (f,s) 

V t , s ip = DZ T -V x cpoZ detD^V*fl>o£= (-3^3,^). V t<8 f. 

Let / = Z(fc) be the non-trivial permutation of {1,2}. Then we write for the fc-th 
component 

detD£ • d W o I = (-1)' • ( 3 f (3 s & • i/,) - 3 s (3 f £, •!/,)). (4.5) 

This equality can easily be derived from the product rule. 

Without loss of generality, we assume that the parametrization of the boundary 
7(s) — £(0,s) is clockwise which implies detD£ > 0. In addition, we have the useful 
relations 

3^=4°^ 3 s ^ = -No^.|3 s f|. (4.6) 

The first equality is clear by the characteristic ODE. The second is a consequence of 
£(f, .) being a clockwise parametrization of the £-level of To- Finally, we remark that 
ip(t,s) is periodic w.r.t. the variable s, since £ is. Thus we have rp(t,a) = ip(t,b). 
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Now we can compute the derivative measure Du, which in both parts is the same 
process. Let cp 6 C] (Cl) be a test function. For the moment we restrict the discussion 
to subsets O a of O which are lower level-sets of To, that is 

n A := n n {x e n : t (x) < a} , 

for < A < 1. Note that C\\ = Cl. When later on we have cp £ C\ (0\E) (respectively 
cp E Cl(Cl\L) 2 ), as is the case for part a), we will choose A big enough such that 
supp cp c Cl\. For part b) we will pass to the limit A — > 1 instead. 

We separately compute the derivatives of U\ and ui- In order to get Dtfi\ we 
have to look at the following integral: 

b A 

J U\(x) d k cp(x) dx — J J vi(t,s) d k cp o c](t,s) detD£(f,s) dt ds 

a 

b A 

(-1)' / / 7*"o(s) ( MUl ■ V) - M*& -rp))dtds. 



n A a 

b A 



a 



By changing the order of integration and using the integration by parts formula for 
functions of one variable, one obtains 



A b 



= (-l) 1 J 7 *u (s) j d t (d s Z r ip)dtds-(-l) 1 j j '>r*uo(s)d e (dtZrip)dsdt 

a « 

b A b 

= (-!)' J 7*u (s) [d s c; r ip]f =0 ds + (-1) 1 J Jd t c; r ipdD 7 *u (s)dt. 



a 



In the last step we used the fact that 7*«o is a periodic BV-function. Because cp has 
compact support in Cl, we have furthermore ip(0,s) = cp o £(0,s) = cp(j(s)) = 0, so 
the result reduces to 



b A 



(-1) 1 J v 1 (\,s) 3 s £ ; (A,s) • t/>(A,s) ds + (-1)' J Jd& ■tpdtdD 7 *u (s) . 



a 



For the vector-valued version we test with cp 6 (O) 2 and obtain 

y div<p(x) dx — j U\(x) dicpi(x) dx+ J U\(x) d2f2(x) dx 

n A n A n A 

b b A 

= -d s cj2(A r s)-ip 1 (X / s)v 1 (A / s)ds - J J -d t c]2-^idtdDy*u (s) 

a a 

fc fc A 

- y 3 s £i(A,s) • ip 2 (A,s) vi(\,s)ds ~ j j 3 f£i ' ip2dtdDy*u (s) 

a a 

J) fe A 

- - y (i/;(A,s),a s ^(A,s)) Pl (A,s)ds - J J (ip^) dtdD<y*u (s) 



a 
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By using the relations (4.6), the last result can be written as 

b 

J m(x) div (p(x)dx = - J (cp,-N) o £(A,s) u x o £(A,s) |3 s £(A,s)| ds 



b A 



/ / (<p,4)o£dfdD 7 * Mo (s). 



a 



Here, the first summand integrates w.r.t. the T^-measure along the A-level of To. For 
the restriction of measures onto A-levels of Tq we will use the abbreviation 

"H 1 LA := H x L{x G n : T (x) = A} . 

In the second integral we rechange variables by pushing forward the product mea- 
sure C 1 ® Dj*Uo with the diffeomorphism g. Let }i denote the pushed-forward mea- 
sure ]i := ^(-C 1 <E> D7*«o)- Then we finally obtain 

J U\(x) div<p(x) dx = — j (cp(x),—N(x)) ui(x) d'H 1 L\(x) 
ci A n 

VW' c o (*)) d F(*)- 



For the derivative of «2 we perform the same steps as above with the integral 

b A 



u 2 (x) d k q>{x) dx = (-1)' J J v 2 (t,s) ( dt(d 8 £i ■ tp) - d e (dt£i ■ ip) ) dt ds . 
n A a 

After changing the order of integration we go on with integration by parts: 

b A A b 

= (-!)'/ Jv 2 (t,s)d t (d s £ r ip)dt ds-(-l) l [ J v 2 (t,s)d s (d t £ r ip) ds dt 

a « 

= (-!)'/ ^2(A,s)a s ^(A,s)-^(A,s)- Jd t vzd s eripdt\ ds 

A & 

- (-1)'/ f-d e v 2 d£ripdsdt 
o « 

f> b A 

= (-1)' jv 2 {Ks) a s ?/(A,s) ■ 0(A,s) ds - I J(-l) l (d s £t d t v 2 - 3 f ft d s v 2 )-ipdtds. 



a 



In the second equality we have used the periodicity of £ w.r.t. s. Because of v 2 
u 2 o £ we can - according to equation (4.5) - substitute the last integrand to get 



& A 



(-1)'/ »2 (A, s) 3 s £;(A,s) -i/>(A,s) ds - J J det D£ ■ \u 2 o £ ■ ip dt ds . 



a 
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By means of the last result and a rechange of variables, we end up with 

J U2{x) div <p(x) dx = — J (q>(x), — N(x)) ii2(x) dH 1 L A(x) 

— J {(p{x) ,V U2(x)) dx , 



when testing with cp £ C^(O) 2 . Adding the partial results for U\ and «2 gives us 
J u(x) div cp{x) dx = — j ((p(x),—N(x)) u(x) dT-L 1 L A(x) 



n 



J (cp{x),c^{x)^ dpi(x) - J (<p(x),V« 2 (x)) 



(4.7) 



dx . 



n, 



Now, we are ready to turn to the proof of the assertions a) and b). 
a) In this part we have Q\E as the domain of u. If we test with q> £ Cj(D\S) 2 , we 
can choose A < 1 large enough that equation (4.7) reduces to 

j u(x) div cp(x) dx = — J (^f(x) , Cq(x) dfi(x) + Vii2(x) dx^j . 
n n 

Consequently, the derivative measure is given by Du = Cq(x) ■ \l + V«2(x) ■ £ 2 - 
What remains to show is the boundedness of || m II bv(0\S) - 

Because uq £ BV (dfi) is a BV function of one variable, uq is essentially bounded 
(see [ , chapter 3.2]). By equation (4.4) then, we obtain a bound on the L 00 -norm 
of u: 



l M llL~(n) < ll M o||i~(3n) + II To II oo ||/o II oo < ||wo||i°°(an) + 



B-m 



For the total variation \ Du\ (Q\£) we estimate both summands of Dm separately, 
beginning with 



(P(x),Cq(x)) d}i{x) 



n 



< 



b 1 



a 



B-m 



d\D 7 *u \(s)- \\f\\ 



>o£,c^-o(f) dtdDj*u (s) 
\Du 



B-m 



m\co ■ 



Clearly, the total variation of this summand is bounded by |cg"(x) ■ }i \ [Cl\L) < 

^ DU ^m^ > wn i cn is the total variation of the boundary data times the bound on 
the arc lengths of characteristics (see lemma 3.1 b) ). 
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The total variation of the second summand is exactly the I^-norm of V u 2 \ 



b 1 



J \Vu 2 {x)\dx = J J |Vm 2 o£ -detD£| dt ds 



n a 

b 1 



a 

b 1 fe 1 

< 

fl a 



1 1 |3 S £| |3*o 2 | dtds + J J \d£\ \d s v 2 \ dt ds 



This step is completed if the last two integrals are bounded. The partial derivative 
of i>2 w.r.t. t is 

3*e 2 =/ o£=^-^o£=/oe.|a f £|. 

Additionally, by using relation (3.8) from corollary 3.3 c), we obtain the bound 

b 1 b 1 

1 1 |3 S £| |3 t z; 2 |dfds = || |3 S £| |3 t £| \fo£\dtds 

a « 

6 1 

< 1^ |detD£|/o£|dfrfs= 1 y |/(x)| rfx < 1 ■ H/IU-r^n) 
fl n 
on the first integral. The partial derivative of c 2 w.r.t. s is 

f 

d,V2 = |(V/ o£(t,s),3 s £(t,s)) rfr. 


With this, we estimate the second integrand by 

1 

|3 f £| |3 s z; 2 | <-^—- I |(V/oo£(t /S ),3 s £(t, S ))| dr. 







Because the latter is independent of t, we obtain 

b 1 b 1 

1 1 13*£| |3 s o 2 |dfds< -g-^-- 1 ||(V/ °£(T,s),a s £(T,s)>| dzds. 

fl fl 

According to the definition of fo in equation (4.2), we expand V/o to 

Moreover, with 3 S £ = N x o £ ■ |8 S £| and |3f£| = 1/ (c, VTo) o £, we end up with 
(V/o o £3,0 = |3 S £||3*£| (( V /^ x )-^- (Voc^c 7 ^!^) o£. 
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Finally, by writing N as N = VT / 1 VT |, we have DN = N L ■ N ±T ■ D 2 T / \ VT 1 
and can express the summand involving D 2 Tq in terms of DN: 

N ^. dn . c = n ±t.D^To. c ^ c t.D^To. n ±_ 
|VT | |VT | 



Hence, there is the integrable upper bound 

iti 

1 



(V/ O o£,a s f)|< ("nv/ll^ + M^dDcl + lDNl) ) of, 



which implies 



/ / |8 f? ||8 s z; 2 |d^ s <^^r 2 (n) + &(||Dc| 



Li(n) + ll DN llLi(n) 



a o 

Combining the partial results we firstly get a bound on the L 1 -norm of V«2 

ll/lioc l|V/JU\ , JJ/I 



and secondly learn that the total variation \Du\ (fi\E) is bounded by Mq\ £. 

b) Now, we want to view u as a BV function on the domain O. Thus, test functions 
stem from C} (O) 2 and we have to study the limit A — » 1 in equation (4.7). In part 
a) we already have bounds on the total variation of the components concerning 
Cq(x) ■ ]i and Vujix) ■ C 2 , which do not depend on A. Hence, these bounds stay 
the same, and we can concentrate on the remainder 



/ 

n 



b 



(cp(x),-N(x)) u(x) dU l L\{x) = J (ip(X,s),d s ^ ± (\,s)^ v(\,s) ds . 



In corollary 3.3 d) we introduced a partition of I, the domain of 7, such that 
f (1, -)\j k+ and f(l, .)|j. are both regular parametrizations of E^, one for the plus 

side and one for the minus side of Ej-, i.e., 

1= U ih+^h-) u /. 

k=l 

The remaining part / comprises all s for which the characteristic f(A, s) hits a 
singular node of E, i.e., a terminal-, branching- or kink node, as A tends to 1. As 
before, let z\, . . . ,z m denote the singular nodes of E. Then we partition / into a 
collection / Zl ,...,/ Z „, by: 

m 

se} z , limf(A,s) =z\ , / = (J / z , . 

A— »1 . 



For an illustration of this decomposition, see figure 4.1. 
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FIG. 4.1. Decomposition of the level line {Tq — A} = £(A, . ) into 4 parts. The left image shoivs the domain 
fi in white, the boundary dCl in blue, and in red the stop set E with terminal nodes Z\ and Z2 ■ The red vector n is the 
normal ofL and defines the plus and the minus side. The black curves are two characteristics, both start at dCl, but 
the upper one approaches JLfrom the plus side while the other comes from the minus side. The right image illustrates 
the decomposition of £(A, . ), the green curve. In black we have characteristics which hit the terminal nodes. For 
each node only the "first" and "last" one of the characteristics hitting the node are painted. In solid blue we have the 
parts 7(/ Z j) and 7(/z 2 ) which comprise all the start points for which characteristics hit terminal nodes of E. The 
complement is painted dashed blue and consists of the two parts 7(/+) and "({]-)■ Here, 7(/+) (j{J-)) contains 
all start points for which characteristics approach JLfrom the plus side (minus side) and end in the relative interior 
o/E. With parameter sets J Zl , ] + , and J_, that correspond to the decomposition of the boundary curve 7, we 
obtain an analogous decomposition of the level line {To = A} = £(A, . ) hf§(X,] Zl ), £(A, J Z2 ), £(A, J + ), £ (A, /_). 
If now A — ► 1, i.e., {Tg = A} — > E, the curves £(!,/+) and £(!,/_) remain and are two parametrizations o/E. 



By these partitions, we decompose the integral 



m n 

<KA,s),a s f x (A,s)) v{X,s) ds = £ / ... ds + £ 
' 1=1/ fc=i 



For the / z ,-summands, we have the estimate 



J (ip(A,8) f d a g x (K8)) v{A,s)ds 




((cp,-N)-u)(x) dU\x) 



?(A,/z,) 



— Il^ll 00 ' ll M llL~(n) ' 



Because the arc f(A, } Zj ) degenerates to the single point z ; (see figure 4.1), the 
right-hand side becomes zero in the limit 

lim^ 1 (?(AJz ( ))=0, 

A— H 

and the contribution of those summands vanishes. 

For the remaining summands we perform only those which go along £(A, + ), 
since the same argumentation applies to the others which go along £(A, ]j. _). For 
£(1, .) : Jk,+ — > £jt the tangent is given by the limit 



d s £(l,s) = lim a s ?(A,s) 
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By requirement 2.4 part 2, the field of normals extends to £/ c . This means 
lim -No £(A,s) = n k og(l,s) s E } k/+ . 



Finally, the second component s(x) of (see corollary 3.3 b) ) extends one- 

sided onto Lfc, an£ l so does u. Let z 6 1L k , with corresponding s+(z) 6 /j. + . Then 
we set 



uUz) := lim o(A,s + (z)) = v(l,s + (z)) . 



Conversely the following relation holds, 

o £(l,s) = lim u o£(A,s) = o(l,s) , s 6 } kA 

A->1 

Now, we can turn to the limit. For abbreviation let 

h(A,s) := ({q>,-N)-u)o{(\,s) . 

Then we have 



J {{cp,-N) -u){x)dU l {x)- j ({q>,n k ) -u+){x) dU l UL k {x) 



(4.8) 



£(*./*,+) 



n 



J h(A,s) ■ |3 s £(A,s)| ds - J h(l,s) ■ |3 s £(l,s)| ds 

Jk,+ Jk,+ 

< I \h(l/s)-h(\,s)\ ■ \d s £(l,s)\ ds 



h,- 



+ Whor'h^n) J \\dsS(l,s)\-\dsZ(K 

Jk,+ 



s)\ \ ds . 



By the extensions of N and u the product \h(l,s) — h(X,s) \ ■ |3 s g(l,s)| tends to 
zero for every s £ J k+ . Furthermore, it has the integrable bound 

1/2(1,5) -h (A, s) | ■ |3 8 £(l,s)| < 2||^or 1 || L ~(n) ' |3 s ?(l,s)| . 

Thus, by dominated convergence, the corresponding integral vanishes as A — >■ 1 . 
By the same argument the second integral tends to zero, too. 
Summarizing, this means 

(q>, -N) ■ u m\x) -»■ J (cp,n k ) ■ u+ dH 1 L2, k (x) , 

(<P,-N) -udH l {x) ->■ / (<p,-n k )- Uk - dH x l-JL k {x), 

as A — »■ 1, and together we obtain the jump part for 

(<p,-N)-ud-H l (x) -> [(4 - Uk -)-(f,n k ) dH l LJL k {x). 



a^Jk,+)u?(KJ k ,-) 



n 
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The bound on the total variation of the jump part is 



J( u k- U k) (<P' n k) dU x UL k {x) 
o 



< 2 ■ ll u !lL~(n) ■ '^ 1 (^/t) / 



whereas ||f>||oo < 1. What remains to show is that the one-sided limits u^,u^ 
defined above are in fact the BV-traces w^T and w^T of u on Ej-. This is true, but 
we postpone this point to lemma 4.2. 
Clearly, the complete additional jump part is given by 

n 

E / ( u k( x ) ~ u k M) (f' n k) ( x ) dT^LZ^x) 



k=1 6 



with 2||«|| L co(q) ■ ~N}(L) as a bound on its total variation. We have shown that u 
is in fact an element of BV(d) with || « Hbv(O) bounded by the given data, and its 



derivative measure reads 

n 



Du = Y,(u+ {x) - u- {x)) n k (x) ■n 1 Li: k + cfr(x)-fi + Vw 2 (x) • £ 2 
k=l 



Theorem 4.1 part a) shows that the candidate solution u belongs to BV(0\E) 
and the derivative Du, on Q\£, is given by Du = Cq(x) ■ u + V 1*2(3:) ■ £ 2 . From 
orthogonality, we infer 

(c (x),Du) = (c (x),X7u 2 {x)) ■ C z , in fl\L 

and as in the classical case, one checks that (co(x), V«2(*)} = fo( x )) the relation be- 
tween backward and forward characteristics according to equation (3.7) is £(Tq(x) — 
t,s(x)) = t](t,x). Now, substitute x = £(t,s) to get £(t — t,s) = rj(t, £(t,s)), differ- 
entiate w.r.t. t, resubstitute and obtain the relation 

-n'(t,x) =D x n{t,x)-c (x) . (4.9) 

The description (4.4) of u 2 , equation (4.9), and the equality (c , VT ) = 1 imply 

Mx) 

(c (i),V« 2 (x)) = fo°v( T o( x )> x ) ~ J (V/ o^(t ( i),);'(t ; i)) dr = f (x) . 

o 

Consequently, the candidate u solves the PDE of problem (2.1). 

In the next lemma we study the traces of the candidate u along level lines of 
To and observe thereby that u satisfies the Dirichlet boundary condition of problem 
(2.1) as BV-trace. 

Lemma 4.2. (Start /restart /stop) 
a) (start): u satisfies the boundary condition w|gn = uq as BV boundary trace, i.e., 



lim -=■ / \ua(z) — u(x) \dx — 
rH>0+ r 2 J 

nns r (z) 

for every Lebesgue point z G 3fi o/uq. 



THE DIRICHLET PROBLEM FOR LINEAR HYPERBOLIC PDES 



21 



b) (restart): for every z 6 Q\£ that corresponds to a Lebesgue point z' of uq, that means 
z' = 7(s(z)) is a Lebesgue point ofu.Q, we have 

lim 1 n [ \u(z) — u(x)\ dx = and lim / \u(z) — u(x)\ dx = . 
r-»o+ r 2 j r^o+ r 2 j 

B>(z) B<(z) 

Here, B<(z) and B> (z) -/or r smaZ/ enough - denote the cut-off disks 
B<(z) := {x G B,(z) : T (x) < T (z)} , B>(z) := {x G B r (z) : T (x) > T (z)} . 



Let T := ^t 0= a « X-level ofT^for some < A < 1 and let T be oriented by N\j. The 
result above means that the traces ui and u^ are identical, thus the restriction u\j is 
well-defined. Moreover, we have u\x G BV(T). 

c) (stop): for every z G £/t t/iflf corresponds to a Lebesgue point z' of uq w.r.t. the plus side 
o/Efc, f/wf means z' = 7(s + (z)) is a Lebesgue point ofuo, we have 



lim \ f \ut(z) — u(x)\ dx = 
r->0+ r 2 J fc 



3r + (z) 



Here, u£(z) is defined by equation (4.8) and B+(z) - for r small enough - is the cut- 
off disk centered at z, restricted to the plus side. Hence, the trace u£ is given by u£. 
Moreover, we have u£ G BV (Sj.). The analogous result holds true w.r.t. the minus side. 

Proof. Let z G fl\E and let (t, o~) be its characteristic coordinates, i.e., z = £(t, cr). 
First, we discuss part b) and in particular the case of Bf" (z). Choose r > suitably 
small. By changing variables we get 

s+W'+(w) 

~2 y |u(z) — m(x)| cix = ~2 J J \v(t, cr) — v(t, s)| detD£ (f, s) dt ds . 

B>(z) s_(r) T 

The transformed integrand can be estimated by 

\v(t,s) - v(t,o-)\ < \-y*u (s) - 7*u (cr)\ + 0[s -cr) + 0(t - T) , 
because for the parts v 2 of v which concern the right-hand side /, we have 



\v 2 {t,s) - v 2 {r,a) 



t i 

j f ot(h,s)dh- Jf og(h,a) dh 



< 



|V/ o£(/z,s*)| |3,£(fe,s*)l \s-<r\ dh 



f o^(h,cr) dh 



0(s-cr) + 0(t-r). 



As the determinant det D£ is continuous and positive, we estimate it by a constant. 
Finally, any level line of Tq is a regular C 1 curve and so the cut-off disk B>(z) tends 
to a half disk, oriented along the tangent 3 s £(t, c) of this curve. This implies the 
asymptotic formulas 

s_(r) = cr- <D(r) s+(r) = a + 0{r) t+(r,s) = T + 0(r) . 
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Armed with the last considerations, we have 

^ J \u{z) - u{x)\ dx < - I \ 1 *u (s)- 1 *u (a)\ds + O(r) . 

B>(z) s_(r) 

The latter expression, as r — » 0, tends to zero for any Lebesgue-point a of 7*«0/ i e-/ 
for any Lebesgue-point z' = 7(c) of «o- 

So far, we have got the assertion for the case of Bf" (z). In the case of Bf (z) the 
one and only difference is that, after having changed variables, we have to integrate 
the £ -variable over the interval [/_ (r, s), t] and to perform the same steps as above. 
Together, this proves the first statement of part b). 

For part a), we have z £ dd and hence fl fl B r (z) = B>(z). So, the argument 
used above shows that u satisfies the boundary condition as B V-trace. 

The first statement of part b) implies that the restriction u\x is well-defined. 
When parametrizing T regularly by £a( s ) := £(^/ s )/ we obtain 

£%u(s) = mo£(A,s) = v(\,s) = 7*m (s) + t> 2 (A,s) s eE. 

Since V2{A, .) is a periodic C 1 function, while J*Uq is a periodic BV function, the sum 
£tw is a periodic BV function, and consequently u | y G BV (T) . This proves the second 
statement of part b). 

Finally, for part c), z belongs to the stop set, z £ £/ c . In this case, we cannot use 
£ because its Jacobian is singular on the stop set. But, as £(1, s) for s £ ]y_ + regularly 
parametrizes Sj; (see corollary 3.3 d)), one can use the local diffeomorphism induced 
by 

y' = -c{y) y(0,s) = f(l,s) , s £ J K+ . 

Then, the same considerations as before yield the assertion. □ 

Remark: part b) of lemma 4.2 is called "restart" because having stopped the 
characteristics at some intermediate A-level T of Tq, the restarted problem 

(c(x),Dzv) =f{x) ■ C 1 in Xt >x\^ , 
w\ T = u\x . 

is of the same type as the initial one. This is, because lemma 4.2 b) ensures that 
u\y £ BV(r). Moreover, when applying the method of characteristics here, then iv 
reproduces u by w = u\ XTg>x . 

5. Uniqueness and Continuous Dependence. We have shown in the previous 
section that problem (2.1) has a solution in BV(d). Its uniqueness can be derived as 
in the C 1 case: by linearity, the solution is unique if the homogenous problem 

(c(x),Du) = , u\ m = 

has a unique solution. With v(t,s) = Mof(f,s) the homogeneous problem in charac- 
teristic variables reads as 



D t v = 0, v(0,s)=0 
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Consequently, v must be constant w.r.t. t and is thus determined by the boundary 
condition. The latter implies that the unique solution of the homogeneous problem 
is zero. 

As the unique solution u, given by equation (4.4), depends linearly on the bound- 
ary data «o and the right-hand side/, the continuous dependence w.r.t. u$ and f is 
already contained in theorem 4.1. What remains to be shown is the continuous de- 
pendence of the solution on the transport field c and the time function T. 

THEOREM 5.1. (Continuous dependence on c and T) Let (c n ) ne t^ be a sequence of 
transport fields and let cbe a transport field. Let {T n )n&i be a sequence of time functions 
and let The a time function. Then, consider the family of linear problems 

(cn(x),Du„) = f(x) ■ jC 2 , u n \ m = u Q , (c n (x),T n (x)) > jS- |VT»(x)| , 
(c(x),Du) = f{x) -C 1 , u\ dn = u , (c(x),T(x))>p-\VT(x)\, 

in n\L. If 

1. the sequence (c n ) ne tj converges uniformly to c: \\c n — c||oo — > 0, 

2. the sequence of derivatives (Dc n ) ne jN is bounded: \\Dc„\\ L i^ < M\ for all neN, 

3. the sequence {Tn)n&i converges uniformly to T: \\T n — T\\oo — >■ 0, 

4. the sequence CVT n ) ne ^ converges uniformly to VT: || VT„ — VT||oo — > 0, 

5. the sequence of derivatives (DN n ) n ^ is bounded: \\DN„\\ L i^ < M-ifor all n G N, 
(N n is the field of normals corresponding to T n ) 

then the sequence of solutions u n converges BV weakly* to u: 

u n — ^ u in BV(Cl) , as n — > oo . 

Proof. To begin with, we show the L 1 convergence of u n to u. Since the solution 
linearly depends on the right-hand side /, it is not restrictive to set f = 0. Again, we 
refer to the scaled PDEs 

(c„ r0 (x), Du n ) = u„ | 3n = "o , 

(c (x),Du) = u\ dn = u , 

the first PDE scaled with 1/ (c B , VT n/ o), the second one scaled with 1/ (c, VTo). 
Thereby T„ q and To denote the transformed versions of T n and T according to equa- 
tion (3.1). 

In the following, t] n denote the backward characteristics associated with the 
transport fields c„ o- Likewise, n denotes the characteristics corresponding to Cq. In 
accordance with equation (4.4), the solutions are 

u n (x) = u (t]„(T n)0 (x),x)) , u(x) = u (n(T (x),x)) . 

We split the proof into seven steps. First step: we consider a lower level set 
O a := {x e O : Tq(x) < A} of T such that 2||«o|| L<»(an) ' ^ 2 (^\^a) < e - Then, we 
can restrict the discussion to Ci\: 

\\ u n ~ W||jrl(n) < \\ u n -M||x,l(n A ) + 2 II "0 1| L~(3H) ' ^ 2 (^\^a) < ||«» ~ "IL^Oa) + £ ' 
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Second step: we assume temporarily that the boundary data is continuous. With 

IK - u llLi(n A ) = / \uo(ijn(T nfl (x),x)) - u (f](T (x),x))\ dx , 
n A 

we will get the assertion by dominated convergence if t] n {T n ,o(x), x) — >■ T](To(x),x) 
pointwise for x 6 Ci\. 

Third step: we choose Uq such that for all n > Uq the stop sets E„ - which cor- 
respond to the time functions T n - are all compactly contained in fi\Cl A which is 
possible because T„ — > T uniformly. Now, we are sure that, for n > uq and x £ Cl\, 
the point x cannot belong to some singular set and thus all curves t] n (.,x) are 
well-defined. 

Fourth step: we need to estimate the difference rj n (.,x) — rj(.,x) of the back- 
ward characteristics. Let x 6 Cl\. The derivative of the difference //„(., x) — rj(.,x) 
obviously satisfies the IVP 

(Vn -y)' = c {V) -CnfliVn) , (in - V) (°' X ) = / 
and hence, integration yields 

f 

(rj n - fj){t,x) = J c (t](r,x)) - c nJ0 {t]„(r,x)) dr . 
o 

The latter equation is valid for t < min{To(x), T n fl(x)}. The first estimate is then 

f 

\tj n - rj\(t,x) < J \c (tl(r,x)) - c (fln(r,x))\ + \c (Tj n (r,x)) - c„ /0 (??„(t,x))| dr . 
o 

According to requirement 2.5 3a) we get a bound |Dc(x)| < M^onO A . Andasimilar 
bound | Dcq (x)\ < M' A holds for the transformed transport field Cq, which implies the 
existence of a Lipschitz constant on Q^. For < t < min{To(x), T n g(x)} all the 
points 7] (t, x) and r\„ (t, x) are located in Q,\ and so the next estimate is 

t 

\rj„ -tj\{t,x)< J L x - \Vn - rj\(r,x) dr + t\\c - c„ /0 || L ~(n A ) / 
o 

and thus Gronwall's lemma (see [ ]) yields 

\Vn ~ rj\(t,x) < t ■ e^ LA -||c - c„ /0 || L co (nA) < C A • \\c - c„ /0 || L ~(n A ) ■ (5-1) 

Because of t < min{To(x), T Hi o(x)} < Tq(x) < A for x G O a , the constant is C,\ = 
A ■ e k . Since we have three times uniform convergence, i.e., T n — > T, VT„ — > VT, 
and c„ — >■ c we obtain uniform convergence of the transformed versions T„ o — »■ 
To, VT H/ o — >■ VTo, and c H/ o c at least on O^. Consequently, we can make the 
difference \t] n — t]\ (f, x) arbitrarily small. 

Fifth step: we estimate the difference \t]n(T n ,o{x),x) — r](To(x),x)\. For abbre- 
viation we set t = Tq(x) and x n = T n fi(x). We want to reuse estimate (5.1), but 
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this estimate is valid only for t < min{T, T n }. So, in the case of T — min{r, r n } we 
proceed as 

\t]n(r n ,x) - rj(r,x)\ < \rj n {r n ,x) - f] n {r,x)\ + \tj„(t,x) - t](t,x)\ 

< II - Cn,0° >7«(-/*)||l~([t,t„]) ■ l T ~ T n\ + C A • ||c - C M/0 1 1 (Q A ) 

- +Ca -ii c °- c "'°ii^(«a)' 

and in the case of T„ — min{T, t„} we supplement the other way round 

\rin{r n ,x) - t](t,x)\ < \rj n (r n ,x) - r]{r n ,x)\ + \rj(r n ,x) - rj{r,x)\ 

< C A ■ 1 1 c - c n/ o||i°°(n A ) + || — c o 77(.,x)||l°°([t„,t]) ■ \t - r„\ 

| T — T M | 

< C A ■ \\c - c„, || L ~ ( n A ) + g. mo ■ 

For each n, the value m n o > denotes the minimum of | VT„ g| according to lemma 
3.1 part a), and because of | VT„ q| — >■ | VTo| locally uniformly we also have m n Q — > 
mo. So, we have everything to conclude \i] n (T n fi(x),x) — t](Tq(x),x)\ — > and back 
to the second step we infer \\u n — »||xi(Q A ) — > 0. 

Sixth step: in the case of general BV boundary data we approximate Ug w.r.t. 
II ■ II i,°°(ao) by a sequence (mq )meN or smooth functions. The difference of the solu- 
tions to 

(c{x),Du m ) =f(x)-C 1 , u m \ dn = u n , (c(x),T(x))>B-\X7T(x)\, 
(c(x),Du)=f(x)-C 2 / u\ m = uo, (c(x),T(x))>B-\VT(x)\, 

is exactly \\u — « m || l co (0) = ll u o — u q lL°°(3n)- The same way, replacing c with c„ and 
T with T„, we obtain ||w„ — M™|| L oo(p) = ||«o — Wo'Hi^sn)- After having fixed an m 
such big that 2||«o — u o ! |lL co (3n) ' ^ 2 (Q) < £, we decompose 



\u - u n \\ L i (n) < 2\\u - Uo "|| L » (an) ■ r 2 (n) + j| u m - iC|| Ll(n) < £ + \\u m - <|| 



and apply the previous steps one up to five to estimate remainder \\u m — "Jf Hi,i(o)- 

Final step: all the functions u, u n are elements of BV (CI). By the uniform bounds 
on the derivatives ||Dc n || L n n j < Mi, HDNnHji/Q} < M2, and the boundedness of 
the sequence ffi n ,o, the sequence of total variations |D«„|(n) is bounded. This is a 
consequence of theorem 4.1, where the bounds ||Dc„|| L i^ < Mi, \\DN n || x, 1 (o) — ^2 
are necessary in the case of non-trivial right-hand sides /. By [2, proposition 3.13], 
the boundedness of the sequence \Du n \ (O)and the L 1 convergence u n — y u together 
imply the BV weak* convergence of u n to u. □ 

The well-posedness of the linear problem (2.1) is complete: we have shown the 
existence of a unique solution in BV (O) and its continuous dependence on all of the 
data. Now, we discuss the quasi-linear problem. 

6. The Quasi-Linear Problem. The subject of this section is to proof the exis- 
tence of a solution to the quasi-linear problem 

(c[u}(x),Du) = f[u](x) ■ £ 2 in Q\L, 

"|an = "0/ (6-1) 
(c[u](x),X7T(x)) >B-\X7T(x)\ , 
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in BV(Q). Here, we allow the dependencies of c and / on the solution u to be of a 
general functional type. 

In the introduction we have pointed out, that we achieve this goal by fixed point 
theory applied to the operator U[v), which denotes the solution operator to the fam- 
ily of linear problems 

(c[v](x),Du) = f[v](x) ■ £ 2 in Q\L , veXcBV(n) 

u\dn = u , (6.2) 

(cM(x),vr(x)> >B-|vr(x)| . 

To be able to pursue this strategy, we have to extend the list of requirements on the 
coefficients c and / by assumptions concerning the functional argument V. 

Requirement 6.1. (Transport fields) Regarding the functional argument, transport 
fields are maps of type 

c : L\n) ->■ C^ClXZ) 2 , with c[ .](x) : L 1 (n) — >• R 2 , 

and we assume them to satisfy: 

1. For fixed v G L 1 (D) the function c[v] : 0\E — > M 2 is a transport field according to 
requirement 2.5. 

2. Uniformity of the unit speed and causality condition: 

a) \c[v](x) \ = 1 for all x G n\E and for all v G 

b) There is a uniform lower bound 6 > such that 

6 < (c[v](x),N(x)) < 1 Vxen\S: and Vi>GL x (n). 

c) Both conditions hold for the one-sided limits ofc[v] on the relatively open C 1 arcs Ejt 
ofZ. 

3. Bounds and continuity: 

a) The map D x c : L 1 (D) — > C(0\£) 2x2 - the derivative of c[v] w.r.t. the variable x - 
is L 1 bounded by 

ll D *cM|| L i { n) < Mi VoeL^O). 

b) c is continuous in the following manner: if v G L 1 (Q) and (o n )nsN * s a sequence in 
L 1 (D) with \\v — i>n|li,i(n) ~~ 0/ f'len the sequence (c[v n ]) n <=f$ converges uniformly 
to c[v], i.e., \\c[v) - c[v n ]\\oo -> 0. 

Requirement 6.2. (Right-hand sides) Regarding the functional argument, right- 
hand sides are maps of type 

f : L J (n) -> C J (n) , with /[•](*) '■ — >■ ]R , 

and zt>e assume them to satisfy: 
a) The map f is bounded by 

\\f[v]\\oz < M 2 Vt?GL 1 (n). 
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b) The map V x f : L 1 (Q) — > C(fi) 2 - the derivative of f[v] zv.r.t. the variable x - is 
bounded by 

llVx/MIU <M 3 VoeL 1 ^). 

c) f is continuous in the following manner: if v G L 1 (fl) and (i'n)nsN zs a sequence in 
L 1 (fl) wzt/i ||p — f»||jLi(o) ~~ ^ ®' the sequence {f[v n ])neN converges uniformly to 
f[v], i.e., \\f[v] -/K]||oo -> 0. 

In addition, we define the subsets of BV(dCl) and BV(fi) with which we will 
work later on. 

Definition 6.3. Let M\, Mz, M3 be the bounds from the requirements stated above. 

a) We denote by 

«B = «8(3n) := {v G BV(9n) : ||o|| L -( 3n ) < M 4 , |Dt>|(30) < M 5 } 

i/ze set of boundary functions. 

b) Let the constants M* 6 1R and M** G IR be given by 

M« := 2 • f M, + -n\JL) + 1*L- + f ^ + #M • £ 2 (H) (6.3) 



M 2 



Mi + IIDNI 



1.1 (n) 



6 3 ■ mo 
We set 

X = X(n) := {v G BV(n) : |M| L i (n) < m* , |Do|(n) < M**} . 

Our aim is to view problem (6.1) as fixed point problem. The next corollary 
justifies the change of our viewpoint. 

Corollary 6.4. For fixed v G L 1 (Cl), problem (6.2) meets the requirements of the 
linear problem (2.1) and its unique solution, which we denote by U[v], defines an operator 
U:L 1 (n)^X. 

Moreover, after restriction to X, this operator is a self-mapping U : X — > X . 

Proof. Let v G L 1 (D) be arbitrary. By requirement 6.1 part 1 the field c[v] is 
a transport field according to requirement 2.5, while f[v] G C 1 (D). Thus, the re- 
quirements of the linear problem (2.1) are satisfied and the linear theory ensures the 
existence of a unique solution U[v] G BV(Q). Hence, the map U : L 1 (D) — > BV(Q) 
is well-defined. 

By theorem 4.1 we get the following estimates for fixed v: 

H^MHi^n) <ll M olli°°(an) + ~ ' 

\DU[v]\(n) < 2 ■ \\U[v]\\ L ^ n) — + [—^~ + -pT^ 

^( n ) + ^*-(ll Dc MllL 1( n) + l|DN|| Ll(n) 
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Plugging into these estimates the uniform bounds Mi, M2, M3 from requirements 
6.1 and 6.2 and the bounds M4, M5 on uq £ 03, it is easy to see that we obtain the 
uniform upper bounds 

||LTfa]|| L i(n) < , |DLTfa]|(n) < M„ . 

Summarizing, the operator U is in fact of type U : L 1 (fl) — > X. 

Because of X C L 1 (fl) we can restrict the domain of U to X, and obtain the self- 
mapping U : X — >■ X. □ 



The final ingredient that we need is the continuity of the operator U. 

Corollary 6.5. (Continuity ofU) The operator U : X — >■ Xfrom corollary 6.4 b) is 
sequentially continuous w.r.t. the BV weak* topology. 

Proof. Let (z>„)„ gN be a sequence in X which tends to v £ X w.r.t. the BV weak* 
topology. Then, we have in particular |]r? — Z7 W llz.'-fO) ~* 0- According to requirement 
6.1, the sequence c„ :— c[v„] (with limit c := c[v]) and the constant sequence T n — T 
of time functions satisfy the requirements of 5.1. So as in section 5, we conclude 
|| U fan] — U[v]\\ — > 0, and the weak* convergence U[v n ] — v Ufa] is a consequence of 
the boundedness of |DUfa] | (fl) < M„. □ 

With the corollaries 6.4 and 6.5, we have sufficient conditions to conclude the 
existence of a solution to problem (6.1). By corollary 6.4, the quasi-linear problem 
is equivalent to the fixed point problem u = U[u]. The set X from definition 6.3, 
is non-empty, convex and, by [ , theorem 3.23] sequentially compact w.r.t. the BV 
weak* topology. Finally, corollary 6.5 yields the sequential continuity of the map 
Lf : X — »■ X. Hence, the existence of a fixed point u = U[u] is the consequence of the 
Schauder-Tychonoff fixed point theorem (see [ , chapter 9.3]). 

So far we have shown the existence of a fixed point, but uniqueness cannot be 
expected as the following example illustrates. 

Example of non-uniqueness. As before let fl\Z be the domain of the PDE and T a 
time function. We consider an almost linear problem, where the transport field does 
not depend on u, with zero boundary condition 

(c(x),Du) =/[w](x)-£ 2 , «| an = 0, (c(x),WT(x)) >p-\VT(x)\ . (6.4) 

Furthermore, we choose a right-hand side of the form 

/fa](*)=s(|fa|| L i ( n)) ' ( 6 - 5 ) 

which is independent of x. Thereby g : ]R — »■ ]R is a continuous bounded function so 
that / satisfies requirements 6.2. 

After having fixed the functional argument of / in equation (6.4), the solution 
U [v] is given by 



T {x) 

U[v](x) = J f [v]on(T,x)dT , 




(6.6) 
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according to equation (4.4). Because the backward characteristic solves t]'(.,x) = 
— Cq o rj ( . , x) and c is normed, we have 

\rj'{.,x)\ = \c ot](.,x)\ = 

And because f[v] does not depend on x, the integrand in equation (6.6) reduces to 

f [v}ot](.,x) =f[v] ■ ^ oyj{.,x) =f[v] ■ \f] (.,x)\ . 

Thus, !i[c](x) = f[v] ■ a(x) — g ( || || x, 1 (O) ) ' ^ as se P ara ted variables and a{x) 
is the arc length of the characteristic t](. ,x) which connects x = rj (0, x) and the point 
i](Tq(x),x) on the boundary dCl. 

If now u is a fixed point of U, it must satisfy 




Consequently, u(x) = a ■ a(x) is a scalar multiple of the arc length function a. Plug- 
ging this into condition (6.7), we obtain by 

K-a(x) =g[K- ||fl||x,i(n)) •«(*) a =g(*- ll«lli,i(o)) 

a scalar fixed point problem for g. Finally, if g is a continuous bounded function with 
many fixed points a, e.g. 

[-1 ,f<-l 
g(t)= It ,-l<f<l , with g(f)=| 
[l ,1<* 

then the operator U has as many fixed points as g. 

7. Conclusion. We have shown in section 4 that the linear hyperbolic Dirichlet 
problem - with interior outflow set Z and mere inflow boundary - admits a solution 
in BV(fi) under the assumptions of section 2. The crucial ingredient here was the 
knowledge of a time function T w.r.t. which the transport field c is causal. Further- 
more, we have shown in section 5 that this solution is unique and depends continu- 
ously on all the data of the problem. 

Certainly, one might ask why we prefer the space BV(Cl), because, if the bound- 
ary data uq were C , one could solve in C 1 (n\S) (this possibility is contained in our 
approach). The advantage of working in BV (O) is three-fold: 

1. Our solution solves the problem in C1\Z, but is in fact an element of BV(Cl) ac- 
cording to theorem 4.1 part b). Therewith we get a description of what happens 
to u on £ in the language of geometric measure theory and can in this sense close 
the gap Z. 

2. The benefit of closing the gap is that the domain X C BV(C1) of the operator U 
from section 6 is compact w.r.t. the BV weak* topology. This fact together with 
the continuous dependence allowed us to conclude the existence of a solution to 
the quasi-linear case by employing the Schauder fixed point theorem. 



(c,VT ) 



f](.,x) 



(c,vr ) 



Ot](.,x) 
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3. A further advantage is that we can allow for more general time functions T and 
stop sets E. For example we can use a X with a saddle node as shown in figure 
7.1. Here, in a first step, one solves from the boundary to the saddle level. The 
saddle node in the middle is the first stop set so to speak. For the second step, we 
consider two decoupled problems on the remaining parts. At the saddle node we 
will in general get a jump discontinuity even if the data on dD, is smooth. That 
means the problems on the remaining parts - (re-)started on the saddle level - are 
supplied with BV-data. For more on this topic see [ , chapter 5]. 

In this paper we have restricted the discussion to two-dimensional domains and 
in the proofs of theorems 4.1 and 5.1 for arguments concerning the boundary data 
we used strongly the fact that BV functions of one variable are essentially bounded. 
In order to generalize to d-dimensional domains, d > 2, with a [d — 1) -dimensional 
boundary, one must require the boundary data to be BV{dd) n L°°(30). 
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